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^ I We present chemical abundance measurements in the Zabs = 3.35045 Damped 

■ Lyman-a (DLA) system observed in the UVES spectrum of the BAL quasar BR 

r~| . 1117—1329. We measure a neutral hydrogen column density N(HI) = 6.9=bl.7x lO^'^ 

Qh| atoms cm~'^ and derive mean abundances relative to solar: [Si/H] = —1.26 ± 0.13, 

O ■ [Fe/H] = -1.51 ± 0.13, [Ni/H] = -1.57 ± 0.13, [Cr/H] = -1.36 ± 0.13, [Zn/H] = 

'. -1.18±0.13, [Al/H] > -1.25, [0/H] > -1.25 and [N/H] < -2.24. This is the third 

^ I measurement of Zn, an element mildly depleted onto dust grain, at z^fes > 3. The 

• iron to zinc and chromium to zinc ratios, [Fe/Zn] = —0.33 it 0.05 and [Cr/Zn] = 

•'— j . —0.18 lb 0.05 demonstrate that the absorber has a low dust content. The nitrogen 

I ratio [N/Si] < —0.98 suggests that the "secondary" N production process is taking 

' place in this DLA. Finally, this absorber does not seem to present a convincing a- 
enhancement as shown by the a over Fe-peak element ratios: [Si/Fe] = 0.25 ±0.06, 
[Si/Cr] = 0.10 ± 0.06 and [Si/Zn] = -0.08 ± 0.06. 
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1 Introduction 



The chemical abundances of the highest H I column density quasar absorbers 
{N{HI) > 2 X 10^° atoms cm~^), the Damped Lyman-a systems (DLAs), 
are commonly expressed using column density weighted measurements. They 
are used as observational tracers of the cosmological evolution of metallicities, 
up to very high redshifts (see for example Pettini et al. 1994, 1997; Lu et al. 
1996; Prochaska & Wolfe 1999, 2000, 2001; Dessauges-Zavadsky et al. 2001; 
Ledoux, Bergeron & Petitjean 2002 and references therein). However, recent 
studies give somehow surprising results: Prochaska & Wolfe (2002) claim no 
evolution of the mean weighted iron metallicity over the redshift range 1.7 — 
3.5, in contradiction with predictions from essentially all chemical evolution 
models. A possible explanation for this result is that iron is not the optimal 
element for tracing metals on cosmological scales, because it is easily depleted 
onto dust grains. Dust depletion greatly complicates the interpretation of the 
metal content of DLAs. Nevertheless at very high-redshift {z > 3.8), there are 
evidences that [Fe/H] is beginning to fall, in such a way that at z = 5 the 
metaUicity is substantially lower than at 2; < 4 (Prochaska, Gawiser & Wolfe 
2001, Songaila & Cowie 2002). 

Pettini et al. (1994) estimated the dust content of DLAs with the help of a 
comprehensive survey of Zn measurements, an element which is known to be 
only slightly depleted onto dust grains, and thus provide an unbiased tracer of 
metallicities. However such measurements are challenging because of both the 
paucity and weakness of Zn features in quasar absorbers. At the moment, only 
two measurements of [Zn/H] have been made in DLAs at z>3 (Molaro et al. 
2000, Levshakov et al. 2001). Here, we present a third detection of Zn in a DLA 
at z~3.35, thus providing clues on the dust- free abundance determination in 
high-redshift Damped Lyman-a systems. We first detail the observational set- 
ups and data reduction processes of the quasar spectrum. We then present the 
analysis and chemical abundance determination for a number of elements in 
the DLA studied, providing a short discussion on the consequences of these 
additional measurements at high-redshift. 



2 Observation and Data Reduction 

BR 1117—1329 is a high-redshift {z ~ 3.96) quasar exhibiting broad absorption 
lines (BAL) discovered by Storrie-Lombardi et al. (1996). The data presented 
in this study were obtained using the Ultraviolet- Visual Echelle Spectrograph 
(UVES) on the 8.2m VLT in April 2000 and April 2002. UVES is particularly 
well suited for the study of Zn at high-redshift thanks to its good sensitivity 
at the longer wavelengths. Three partially overlapping setups (480-1-710-1-725) 
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Fig. 1. Voigt profile fits to the Lyman series lines (Ly-a, Ly-/? and Ly-7) of the 
z = 3.35045 DLA. The vertical bar in each panel indicates the wavelength cen- 
troid of the component used for the best fit, which is shown as a solid line to- 
gether with 1 a errors (dashed lines). The resulting total HI column density is 
iV(HI)= 6.92 ± 1.7 X 10^0 atoms cm^^. 

were used to provide a complete coverage of the quasar spectrum ranging from 
~ 4195 A to 9200 A. 

The individual spectra were reduced using the UVES data reduction pipeline 
of the ESO MIDAS package (see Ballester et al. 2000 for a detailed descrip- 
tion). The steps of the data processing include wavelength calibration, order 
extraction and flat-fielding. Once extracted, the individual frames were cor- 
rected to a vacuum heliocentric scale and combined together resulting in a 
spectrum with signal-to-noise ratio of about 40 per pixel at ~ 7000A. The 
quasar continuum was then fitted to the quasar spectrum using a spline func- 
tion. 
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Fig. 2. Comparison between the spectrum of BR 1117—1329 (solid) and that of 
another high-redshift quasar taken as a reference (dashed). The superposition of the 

two spectra allows to check for possible atmospheric contamination in the observed 
wavelength range 8968A-8993A. The Zn II A 2062 and Cr II A 2066 arc contaminated 
by sky features while the Cr II A 2062 line appears clear from any contamination. 

A careful inspection of the spectrum bluewards of the Lyman-a emission 
reveals a DLA absorption in a region free of BAL features. The redshift 
(zabs — 3.35045) is well defined by the presence of higher Lyman series lines 
and heavy element absorption lines including some from Zn ll. 



3 Analysis 



3.1 Column Densities Measurements 



The H I column density of this system is measured by fitting a theoretical Voigt 
profile to the Lyman series. The quasar continuum was fitted with a spline 
function in MIDAS. As always in DLAs at high-redshift, the positioning of the 
continuum over the strong absorption lines is the main source of uncertainties 
in the H I column density measurement. The use of several lines to perform the 
fit gives confidence in the result however. The fitting of the profile is performed 
with VPFIT ^ (Webb et al. 1987) and the resulting fit is presented in Figure 
1 for the first three lines of the Lyman series. The total H I column density is 
A^(H l) = 6.92 ± 1.7 X 10^0 atoms cm-^. 

Wc derive the column density of the ions associated with the DLA by fitting 



3 VPFIT was written by R.F.Carsweh, J.K. Webb, M.J. Irwin, A.J. 
Cooke. More information about VPFIT are available at following URL: 



http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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Relative velocity (km/s) 

Fig. 3. Velocity profiles (km s~^) of the metals in the DLA absorption system with 
the fit overplotted (light colored solid line) . The spectra are normalized to unity. The 
zero velocity is fixed at z^fes = 3.35045 and the dashed lines indicate the positions 
of the 3 components. 
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Table 1 

Summary table of the oscillator strength values used in this study. 



Ion 




fosc 


References 


Si II A 1526 


0, 


.11000 


Spitzer & Fitzpatrick (1993) 


Si II A 1808 


0, 


.00218 


Bergeson & Lawler (1993b) 


Fe II A 1608 


0, 


.05800 


Bergeson et al. (1996b) 


Fe II A 1611 


0, 


.00130 


Welty et al. (1999) 


Ni II A 1454 


0, 


.03230 


Fedchak et al. (2000) 


Ni II A 1709 


0, 


.03240 


Fedchak et al. (2000) 


Cr II A 2062 


0, 


.07796 


Bergeson & Lawler (1993a) 


Zn II A 2026 


0, 


.48860 


Bergeson & Lawler (1993a) 


Al II A 1670 


1, 


.83300 


Morton (1991) 


I A 971 


0, 


.01214 


Morton (1991) 


I A 1039 


0, 


.00920 


Morton (1991) 


N I A 1200.2 


0, 


.08849 


Morton (1991) 


C IV A 1548 


0, 


.19080 


Morton (1991) 


C IV A 1550 


0, 


.09522 


Morton (1991) 


Al III A 1854 


0, 


.56020 


Morton (1991) 


Al III A 1862 


0, 


.27890 


Morton (1991) 



theoretical Voigt profiles to the observed features minimizing x^- For this 
purpose we make use of the FITLYMAN package in the MIDAS data reduction 
software. Only the lines free of telluric contamination and blending are used 
to determine the column densities of the elements. In particular, we took great 
care to check possible atmospheric contamination in the Cr ll A 2062 region 
(see Figure 2). 



3.2 Results 



The column densities obtained from the Voigt profile fits are used to derive 
abundances of Fe, Si, Ni, Zn, Cr and to place a lower limit on the abundance 
of Al and O and an upper limit on the abundance of N. All the low ionization 
species were fitted together with 3 components (sec Table 2). We notice that 
all the lines can not be fitted with a single Doppler parameter. However a 
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Tabic 2 

Results of Voigt profile fits to low ionization species in the z = 3.35045 DLA system. 



Component 




Ion 


[km/s] 


log N(X)'' 
[cm~^] 


A 


3.350265 


Fell 1608-1611 


6.5±0.5 


14.39±0.01 






Sill 1526-1808 


4.0±0.4 


14.34±0.05 






Nill 1454-1709 


6.5±0.5 


13.07±0.02 






CrII 2062 


6.5±0.5 


12.76±0.02 






Znll 2026 


6.5±0.5 


11.79±0.03 






AlII 1670 




>13.85^ 






01 971-1039 




>15.85'= 






NI 1200.2 




<13.75^ 


B 


3.350452 


Fell 1608-1611 


4.9±0.7 


14.50±0.03 






Sill 1526-1808 


14.8±0.1 


15.04±0.05 






Nill 1454-1709 


4.9±0.7 


13.14±0.02 






CrII 2062 


4.9±0.7 


12.80±0.02 






Znll 2026 


4.9±0.7 


11.86±0.02 






AlII 1670 




>13.65'= 






01 971-1039 




>16.15'= 






NI 1200.2 




<14.39^ 


C 


3.350650 


Fell 1608-1611 


7.7±0.5 


14.09±0.01 






Sill 1526-1808 


15.0±0.1 


13.51±0.05 






Nill 1454-1709 


7.7±0.5 


12.88±0.04 






CrII 2062 


7.7±0.5 


12.30±0.05 






Znll 2026 


7.7±0.5 


11.67±0.04 






AlII 1670 




>11.80'= 






01 971-1039 




>14.50'= 






NI 1200.2 




< 13.55'- 



" The good fit is achieved by assuming two different Doppler values respec- 
tively for Fe, Ni, Zn, Cr and Si, Al, O, N. The latter set of lines is broader 
than the other one. This difference is possibly due to the presence of extra 
components in these lines. 

^ Limits are 3-a. Using same b-parameter as in Si ll. 
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good fit is achieved by assuming two different values respectively for Fe, Ni, 
Zn, Cr and Si, Al, O, N. The latter set of fines are sfightly broader than the 
other one. This difference is possibly due to the presence of extra components. 
It must be noted that this does not affect the column density determination 
as most of the species have at least one transition which is optically thin. An 
exception may be Si ll for which the weakest observed transition Si ll A 1808 is 
partially saturated (although not going to zero). Errors indicated in the tables 
are statistical errors in the fit. However, there are also systematic errors not 
accounted for here but which are smaller than the statistical errors. 

The metal content with respect to solar values is expressed in the usual 
manner: [X/H] = \og[N{X) /N{H)]dla - \og[N{X) /N{H)]q, assuming that 
N{X) = N{XII) and N{H) = N{HI), where the solar values are taken 
from Grevesse & Savage (1998) and Holweger (2001; see Table 4). Silicon 
abundance is derived from the slightly saturated (but not going to zero) 
Si II A 1808 transition and the saturated Si ll A 1526 transition. This re- 
sults in a [Si/H] value of —1.26 ± 0.13. The error in the Si ll column den- 
sity also includes the uncertainty in b's determination. Similarly the iron 
abundance is obtained from the optically thin Fe II A 1611 and saturated 
Fe II A 1608 transitions. We derive [Fe/H] = -1.51 ± 0.13. Two Ni II lines 
(Ni II AA 1454, 1709) are available for abundance determination. We measure 
the column density: [Ni/H] — — 1.57 ± 0.13. The Cr ll fines are detected 
(Cr II AAA 2056,2062,2066), but only one line is suitable for abundance de- 
termination (see Figure 2). The abundance derived from the Cr IlA 2062 is 
[Cr/H] = —1.36 ± 0.13. After checking for possible atmospheric contami- 
nations at these wavelengths, we find the Zn ii A 2026 line allows for the 
determination of the Zn content of this DLA. The new measurement of Zn at 
Zabs > 3, [Zn/ H] = —1.18 ± 0.13, is higher than other measurements at these 
redshifts. The Al IlA 1670 line is heavily saturated and hence only a lower 
limit can be derived from this absorption feature: [Al/H] > —1.25. Similarly, 
several O I lines are detected but all are saturated. Using O 1 A 971 and 01 
A 1039, we derive: [0/H] > -1.25. Finally, the two red most triplets of N I 
are observed but are severely blended. Only one of the lines at our disposal is 
suitable to determine an upper limit to the N I abundance in this system: N 
I A 1200.2. We derive: [N/H] < -2.24. 

We also fit the high-ionization species in the DLA. The resulting column 
density estimates are summarized in Table 3 together with the redshift and 
Doppler parameters associated with each individual component. The C IV 
doublet (C lvAA1548,1550) is spread over ~ 300 km s^^ It is well fitted by 
5 components (Figure 4). We note that there is a velocity shift between the 
centroid of the low-ionization fine profiles (Si ll A 1808, Fe ll A 1608, etc) 
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Fig. 4. Velocity profiles (km s^^) of heavy element lines in the DLA absorption 
system with the fit overplotted (light colored solid curve). The spectra are normal- 
ized to unity. The zero velocity is fixed at the redshift of the absorption feature: 
zabs = 3.35045. 



and the one from the C iv doublet. This is illustrated in Figure 5 where the 
Fe 11AI6O8 line profile is plotted together with the two members of the C IV 
doublet on the same velocity scale. 

Finally, the Al ill doublet (Al lllAA1854,1862) is fitted with two components 
which have their positions and 6-values fixed to those of Si ll. Indeed, the 
velocity structure of Al ill is known to follow closely the one from singly 
ionized species in DLAs (Lu et al. 1996, Prochaska & Wolfe 1999, Prochaska 
& Wolfe 2000). This observational fact is not straight forward to explain: the 
ionization potential of Al ll being greater than the one of hydrogen, Al ill is 
not expected to be present in the regions dominated by neutral hydrogen (i.e. 
DLA) although the atomic physics of aluminum species is not well known (see 
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Tabic 3 

Results of Voigt profile fits to high-ionization species in the z = 3.35045 DLA system. 



Component 




Vel 


Ion 


bthermal 

[km/s] 


log N(X) 
[cm~^] 


A 


3.350745 


-25 


CIV 


20.9+0.4 


13.73+0.01 


B 


3.351459 


+24 


CIV 


16.6+0.9 


13.14+0.02 


C 


3.352319 


+83 


CIV 


12.4+1.0 


12.95+0.02 


D 


3.353086 


+136 


CIV 


6.0+1.2 


12.78+0.07 


E 


3.353214 


+145 


CIV 


24.5+0.4 


13.73+0.01 


a 


3.351072 


-6 


Aim 


8.5+2.4 


12.17+0.02 


b 


3.351246 


+7 


Aim 


5.2+1.4 


11.92+0.03 



Petitjean et al. 1994). Several studies have tried to explain the observations 
(Howk & Sembach 1999, Izotov, Schaerer & Charbonnel 2001, Vladilo et al. 
2001). 



4 Discussion and Conclusion 

In the DLA studied here, we observe a velocity shift between the low and 
the high ionization ionic features. The fact that the high-ionization profiles 
show much more disturbed velocity structure than the low-ionization profiles 
has already been noticed in the past (Lu ct al. 1996; Lcdoux et al. 1998). 
Haehnelt, Steinmetz & Ranch (1998) have suggested that such feature could 
be the signature of merging protogalactic clumps. Indeed they use numerical 
simulations to model the Si 11AI8O8 and C lvA1548 line profiles and note that 
the absorption features vary independently in the high ionization and low 
ionization species since the C IV absorption arises mainly from the warmer 
gas outside the self-shielding region of DLAs. Nevertheless, other numerical 
simulations fail to reproduce the DLA kinematics in a self-consistent manner 
(e.g. Prochaska & Wolfe 2001) and actually require that a significant fraction 
of DLAs have Vdrc ~ 150 km s~^ (Mailer et al. 2001) even within the Cold 
Dark Matter cosmology. In other words, if one wants to interpret the global 
kinematics of the whole DLA population with a simple model, the kinematics 
of the strongest systems might rule out a global dwarf galaxy scenario if one 
adopts Vcirc < 100 km s~^ as the dwarf criterion. Nevertheless, note that in the 
case of BR 1117—1329, the neutral component has a velocity width of about 
50 km s^^ and that the strongest C IV components are seen on both side (at 
—50 and +130 km s^^) of the low-ionization profile. This is suggestive of a 
galactic wind from a dwarf galaxy. 
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Fig. 5. Fe II A 1608 and the C IV doublet (A 1548 and A 1550A) overplotted on the 
same velocity scale. This figure illustrates the shift between the center of the low 
and the high ionization ionic features in the DLA. The zero velocity is fixed at the 
redshift of the absorption feature: Zabs — 

3.35045. 

The detection of N I in this DLA system allows us to probe the chemical history 
of the system. Indeed, the [N/0] ratio is useful to disentangle the contribution 
from the "primary" and "secondary" production processes for N enrichment. 
The hmits obtained here suggests that [N/0]< —0.9 which is shghtly lower 
than what one expects (i.e., ~ —0.6) if the primary process dominates (Vila- 
Costas & Edmunds, 1993). Like in most of the DLAs the [0/N] limits for the 
measured a— element metallicity is consistent with the delayed production of 
N compared to O and inconsistent with the primary production of N from the 
massive stars in the low metallicity gas (Lipman, Pettini & Hunstead, 1995; 
Lu et al. 1998; Centurion et al. 1998; Pettini et al. 2002 in press; Prochaska 
et al. in press) . 

The abundances of the DLA have been studied in detail. In this DLA, we find 
a zinc abundance shghtly higher than the iron one: [Fe/Zn] — —0.33 ± 0.05 
which suggests that the amount of dust present in this absorber is rather 
low. This is further supported by the chromium over zinc ratio [Cr/Zn] = 
—0.18 ± 0.05, another tracer of dust. Thus in this DLA the abundances of 
refractory elements (such as Fe and Cr) are in line with the one from the non- 
refractory element Zn. In addition, Al lll/Al ll is less than 5%. This means 
that the effect of ionization correction in the measurement of metallicity in this 
system is almost negligible (Vladilo et al. 2001) and the absolute abundances 
derived here do not suffer severe biases due to dust depletion. 

We derive various a-chain element to Fe-peak element ratios: [Si/Fe] ~ 
0.25 ± 0.06, [Si/Cr] = 0.10 ± 0.06 and [Si/Zn] = -0.08 ± 0.06. Thus the 
abundance pattern in this system is consistent with no a— enhancement with 
very little depletion in Fe and Cr with respect to Zn. The results indicate that 
the Type II SNe enrichment process does not dominate in this system, pos- 
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Table 4 

Summary table of the column densities of the low and high ionization species in the 
Zabs = 3.35045 DLA. 



Ion 


log [N(X)/N(H)]0 (Ref) 


log N(X) 


[X/H] 


H I 




20.84±0.12 




Fe II 


-4.50 (1) 


14.83±0.03 


-1.51±0.13 


Si II 


-4.45 (1) 


15.13±0.05 


-1.26±0.13 


Ni II 


-5.75 (1) 


13.52±0.03 


-1.57±0.13 


Cr II 


-6.33 (1) 


13.15±0.03 


-1.36±0.13 


Zn II 


-7.40 (1) 


12.26ib0.03 


-1.18±0.13 


Al II 


-5.53 (1) 


>14.06 


> -1.25 


01 


-3.26 (2) 


>16.33 


> -1.25 


NI 


-4.07 (2) 


<14.53 


< -2.24 


C IV 




14.13±0.01 




Al III 




12.36±0.03 





References: 
[1] Grevesse Sz Savage 1998. 
[2] Holweger 2001. 

sibly because the Type I SNe enrichment process is already important. The 
fact that we observe the DLA when Type I SNe processes start to dominate 
over Type II SNe can not easily be interpreted since the turn-over depends 
on both the star formation rate and the initial mass function (Matteucci & 
Recchi 2001). 

We derive a Zn abundance, [Zn/ H] — —1.18 ± 0.13 at Zabs = 3.35, which is 
higher than the very few measurements available at these redshifts (corrected 
to the same solar value as use in the present study): Molaro ct al. (2000) de- 
rives [Zn/H] = -2.00 ±0.1 in a DLA at Zabs = 3.39 towards QSO 0000-2620 
and Levshakov et al. (2002) find [Zn/H] = -1.43 ± 0.08 in an absorber at 
Zabs — 3.02 towards QSO 0347—3819 although the mctallicity in this absorber 
is uncertain and could be larger (sec Lcdoux, Srianand & Pctitjean 2002). The 
resulting column density weighted mean at 2; > 3 is [< Zn/H^LA >] = —1-63, 
i.e. dominated by the highest HI column density, which here corresponds to 
the lowest Zn measurement from Molaro et al. (2000). For comparison, Pet- 
tini et al. (1997) find [< Zn/HoLA >] = -1-18 from z = 0.69 to ^ = 3.39 
(where the measures at 2; > 3 are upper limits). Although statistics arc small, 
our result suggests an increase of metallicities from 2; > 3 to present times. 
Previous compilations of Zn measurements at all redshifts have shown that 
the column density weighted metallicity does not evolve with time (Pettini et 
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al. 1997, Vladilo et al. 2000, Prochaska k Wolfe 2002). These authors argue 
however that this cannot be directly interpreted as a lack of evolution of the 
cosmic metallicity of DLAs since measurements are still affected by low num- 
ber statistics and possible bias against high column density, high metallicity 
absorbers (as first pointed out by Fall & Pei 1993 and Boisse et al. 1998). 
The presence of this bias is controversial at the moment (Ellison et al. 2001, 
Prochaska & Wolfe 2002, Petitjean et al. 2002) and many more measurements 
of undepleted metals in quasar absorbers are needed at 2; > 3 before any 
conclusion can be drawn. 
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